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2001-2002  Annual  Progress  Report  (Year  2) 

This  report  presents  the  specific  aims  and  accomplishments  of  our  breast  cancer  research 
project  during  the  second  year  of  funding  sponsored  by  the  U.S.  Department  of  die  Army.  It 
covers  our  activities  from  September  1 ,  200 1  to  October  3 1 , 2002. 

Introduction 

The  overall  goal  of  this  research  project  is  to  develop  and  evaluate  a  new  approach  to 
monitoring  of  oxygenated  hemoglobin  concentration  (Hb02)  of  breast  tumors  under  respiratory 
interventions  using  near  infrared  (NIR)  spectroscopy  and  imaging  techniques.  Furthermore,  we 
wish  to  compare  and  validate  the  optical  method  with  concurrent  measurements  of  tumor  oxygen 
tension  with  the  use  of  19F  EPI  relaxation  of  hexafluorobenzene  (HFB).  Once  we  have  rigorously 
established  the  integrated  approach  to  tumor  physiology,  we  will  explore  the  influence  of  three 
interventions  expected  to  modify  tumor  physiology.  A  better  understanding  of  the  interplay  of 
these  parameters  in  the  natural  history  of  a  tumor  could  lead  to  enhanced  therapeutic  approaches 
and  provide  a  novel  diagnostic/prognostic  tool  for  breast  cancer  research  and  clinical  practice. 

The  overall  project  has  four  specific  aims: 

Aim  1:  To  evaluate  a  single-channel,  dual  wavelength,  NIR,  frequency-domain  oximeter  and  the 
algorithms  for  obtaining  tumor  Hb02  against  tumor  p02  measured  by  ,9F  EPI  relaxation  of  HFB. 
Aim  2:  To  modify  the  single-channel  system  into  a  3-channel  NIR  system. 

Aim  3:  To  investigate  heterogeneity  of  Hb02  in  breast  tumors  using  the  3-channel  NIR  system. 
Aim  4:  To  study  the  influence  of  three  interventions  on  Hb02  and  p02  of  the  tumor 
simultaneously. 

Specifically,  Task  2  was  planned  for  months  12-18  to  accomplish  Aim  2,  and  Task  3  was 
planned  for  months  18-24  to  accomplish  Aim  3: 

Task  2:  To  modify  the  single-channel  system  into  a  3-channel  NIR  system  so  that  three 
measurements  in  both  transmittance  and  reflectance  geometry  can  be  taken  simultaneously, 
(months  12-18); 

Task  3:  To  investigate  heterogeneity  of  Hb02  in  the  tumor  vasculature  of  the  rat  breast  tumors 
using  the  3 -channel  NIR  system  (months  18-24). 

Body  of  the  Report 

The  PI  and  her  graduate  students  have  made  significant  efforts  to  accomplish  the  tasks. 
We  have  obtained  an  8-channel  NIR  DC  (i.e.,  not  amplitude-modulated)  system,  which  allows 
multi-channel  simultaneous  measurements  in  different  locations.  We  have  designed  and 
developed  a  breast  tumor  dynamic  phantom  in  order  to  perform  3 -channel  experiments  with  the 
multi-channel  NIR  system.  We  have  conducted  computer  simulations  to  further  support  our 
hemodynamic  phantom  measurement  and  to  understand  the  effects  of  tumor  heterogeneity  on  the 
multi-channel  NIR  readings.  Furthermore,  we  have  implemented  coupling  components  in 
preparation  for  the  simultaneous  measurements  of  NIR  and  19F  EPI  relaxation  of  HFB. 
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1.  Tumor  Hemodynamic  Phantoms 


1.1  Motivation  and  the  Phantom  Development 

Our  previous  in  vivo  animal  study  has  clearly  demonstrated  that  cabogen  (95%  02  +  5% 
C02)  inhalation  can  improve  the  vascular  oxygen  level  of  breast  tumors  monitored  by  NIR 
spectroscopy  (NIRS)1.  The  observed  changes  in  tumor  vascular  oxygenation  exhibited  a  rapid 
increase,  followed  by  a  gradual  but  significant  increase,  in  response  to  carbogren  intervention. 
We  have  earlier  developed1  a  mathematical  model  to  interpret  this  bi-phasic  tumor  hemodynamic 
feature:  namely,  tumor  vasculature  is  comprised  with  a  well-perfused  and  poorly  perfused  region. 
To  provide  actual  support  for  this  model  and  to  further  study  heterogeneity  of  the  breast  tumor 
vasculature,  we  designed  and  developed  a  breast  tumor  dynamic  phantom  in  order  to  perform  3- 
channel  experiments  with  the  multi-channel  NIRS  system. 

Gel  phantom  was  fabricated  by  adding  50  g  of  gelatin  powder  (Sigma,  Gelatin  Type  A, 
St,  Louis,  MO)  into  a  boiling  intralipid  solution,  which  consisted  of  350  ml  of  water  and  200  ml 
of  intralipid  (Intralipid®  20%,  Baxter  Healthcare  Corp.,  Deerfield,  IL  20%).  The  use  of  intralipid 
is  to  simulate  light  scattering  in  tumor  tissues.  The  solution  was  stirred  thoroughly  until  the 
gelatin  powder  was  dissolved  in  the  intralipid  solution  completely,  and  then  the  solution  was 
cooled  down  at  room  temperature.  Just  before  the  solution  started  to  be  solidified,  it  was  poured 
into  a  cylindrical  container  (diameter  =  4,5  cm  and  height  =  3.5  cm),  which  contained  two  small 
tubing  networks  used  to  simulate  vascular  networks.  After  completely  being  cooled  down,  the 
solution  became  a  soft  gelatin  phantom  with  the  two  networks  imbedded  inside. 

The  two  imbedded  tube  networks  were  used  to  simulate  two  distinct  perfusion  regions  in 
a  breast  tumor.  We  wound  a  small-diameter  (I.D.=0.52  mm)  tube  and  an  intermediate-diameter 
tube  (I.D.=0.86  mm)  to  simulate  well-perfused  and  poorly  perfused  region,  respectively,  on  a  big 
diameter  core  tube  (O.D,=14  mm).  An  example  of  such  a  tubing  network  is  shown  in  Figure  1, 
A  stream  of  ink  solution  was  used  to  go  through  the  network  to  simulate  a  blood  flow  through  the 
tumor  vasculature.  In  this  way,  we  could  mimic  a  hemodynamic  process  of  a  breast  tumor  by 
changing  the  size  or  length  of  the  small  tube  or  by  controlling  the  flow  rate  of  the  ink  solution. 


Figure  1.  Schematic  diagram  for 
one  imbedded  tubing-network:  the 
small  tube-network  wound  outside 
simulates  blood  vessels  within  a 
breast  tumor. 
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1.2  Experimental  Setup  for  the  Phantom  Study 

Experimental  setup  for  the  tumor  phantom  study  with  the  multi-channel  NIR  system  is 
shown  in  Figure  2.  As  shown  here,  three  NIR  optical  detectors  were  placed  on  the  side  of  the 
cylindrical  phantom;  two  were  placed  in  reflectance  geometry  with  respect  to  the  light  source, 
whereas  the  other  was  in  the  transmission  mode.  With  such  a  setup,  the  two  detectors  in  the 
reflectance  mode  (D2  and  D4)  would  get  the  NIR  signals  from  only  one  of  the  individual  tubing- 
networks,  and  the  third  detector  (Dl)  in  the  transmission  mode  would  detect  the  signal  passing 
through  both  the  networks. 

During  the  measurement,  the  network  tubes  were  filled  with  water,  and  a  diluted  black 
ink  solution  was  used  to  simulate  the  light  absorption  effect  of  blood  and  was  injected  into  only 
one  tube,  followed  by  a  washout  with  water.  Then,  the  experiment  was  conducted  by  injecting 
the  ink  solution  into  the  second  tube,  and  also  followed  by  a  washout  of  water.  Finally,  the  ink 
solution  was  injected  into  both  of  the  two  tubes  at  the  same  time  with  the  following  washout 
water.  Various  flow  rates  (from  2  ml/hr  to  40  ml/hr)  were  used  to  simulate  different  perfusion 
rates  in  tumors.  Changes  of  light  intensities  were  recorded  with  the  computer  through  the  entire 
experiment,  and  the  time  constants  for  the  dynamic  changes  in  light  absorption  were  calculated 
by  fitting  the  data  using  Kaleidagraph  software  (Synergy  Software,  Reading,  PA), 
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Figure  2.  Experimental  setup  for  the  dynamic  phantom  study 

As  shown  in  Figure  2,  detector  2  (D2)  was  placed  on  the  side  near  the  large  diameter 
tube,  and  detector  4  (D4)  was  put  near  the  small  diameter  tube.  Detector  1  (Dl)  was  placed  just 
across  the  source  about  at  the  center  between  the  two  hidden  networks.  Therefore,  we  expected 
that  D2  and  D4  would  get  signals  from  only  large  diameter  tube  and  small  diameter  tube, 
respectively,  and  Dl  would  obtain  signals  from  both  of  the  networks  since  it  had  similar 
distances  to  both  of  them. 
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1,3  Experimental  Results  of  the  Phantom  Measurements 

Figure  3  shows  the  absorption  changes  (Aga,  wavelength=730  nm)  of  tumor  phantom  at 
three  different  positions  when  the  ink  solution  was  injected  into  the  network  tubes.  The  first  step 
for  this  experiment  was  to  inject  the  ink  solution  into  only  the  large  diameter  tube  with  a  flow 
rate  of  20  ml/hr.  As  we  can  see  from  figure  3,  the  reading  from  D2  (pink  trace)  showed  a  largest 
increase  in  change  of  absorption,  and  D4  (green  trace)  showed  the  smallest  increase  in  Aga.  The 
Aga  increase  measured  from  D1  (blue  trace)  was  between  those  obtained  from  D2  and  D4.  A 
similar  pattern  is  also  shown  at  step  2,  where  the  ink  solution  was  injected  only  into  the  small 
diameter  tube  with  the  same  flow  rate  (20  ml/hr).  In  this  case.  Figure  3  demonstrates  that  the 
signal  from  D4  had  a  largest  increase  in  Aga,  and  D2  detected  the  smallest  Aga  increase.  This 
clearly  demonstrates  that  signal  from  heterogeneity  can  be  detected  stronger  if  it  is  near  the 
detector. 


Figure  3.  3-channel  NIRS  results  measured  on  a  tumor  dynamic  phantom  with 
two-size  tubing  networks.  Three  traces  represent  the  readings  at  D1  (in 
transmission  mode  and  located  between  the  two  tube  networks),  D2  (near  the 
large  tube),  and  D4  (near  the  small  tube). 


The  third  step  in  the  measurement  was  to  inject  the  ink  solution  into  both  of  the  tubes 
simultaneously  with  the  same  flow  rate  as  before  (20  ml/hr).  At  this  step,  we  can  observe  two 
features:  1)  the  time  profile  taken  at  D4  (green  trace)  has  a  faster  transition  time  than  that  at  D2 
(pink  trace),  and  2)  the  time  profile  taken  at  D1  (blue  trace)  has  a  bi-exponential  characteristic, 
which  is  very  similar  to  those  we  often  observed  in  our  breast  tumor  dynamic  measurements1. 
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To  understand  these  phantom  measurements,  let  us  consider  the  velocities  at  each  tube. 
Even  though  the  flow  rate  was  the  same  for  both  of  the  tubes,  the  velocities  of  the  ink  solutions 
in  the  two  tube  networks  were  different  due  to  the  different  diameters.  The  solution  in  the  small 
diameter  tube  had  a  velocity  2.74  times  faster  than  that  of  the  large  diameter  tube,  based  on  the 
following  relationship: 

V small  _  r/arge  _  0.43 2  _  2  -74 

^  arg  e  rsmal^  0.262 

where  Vsmau,  and  Viarge  are  velocities  of  the  ink  solution  in  the  small  and  large  tubes,  respectively, 
and  rsman, and  riarge  represent  the  radii  of  the  small  and  large  tubes,  respectively.  This  velocity 
difference  may  be  the  reason  why  the  transition  time  is  faster  in  the  small  tube  than  that  in  the 
large  tube.  Furthermore,  since  D1  was  in  the  transmission  geometry  and  was  located  near  both  of 
the  tube  networks,  the  signal  obtained  at  D1  should  reflect  the  dynamic  changes  occurred  in  both 
of  the  tubes.  We  expect  that  the  bi-exponential  feature  recorded  by  D1  results  from  a 
superposition  between  the  two  different  dynamic  transitions  at  the  two  networks. 

The  phantom  data  presented  here  basically  demonstrate  that  we  can  experimentally  mimic 
the  bi-exponential  behavior  with  two  different  velocities,  or  perfusion  rates,  from  two  different 
diameter  networks.  This  knowledge  may  help  us  understand  the  bi-phasic  features  observed  in 
the  breast  tumor  vascular  oxygenation  measurements  under  carbogen  inhalation.  Furthermore, 
the  next  example  demonstrates  that  the  bi-exponential  dynamics  may  also  result  from  different 
flow  rates,  given  the  same  diameters  for  the  two  tube  networks. 


Step  1  Step  2  Step  3  Step  4 


Top  Tube  I  bottom  Tube  I  Both  Tubes  |  t  I  Both  Tubes 


no  ml/hr)  I  (20  ml/hr) !  |  (20  ml/hrt  I  !  I  (5  +  20ml/hri 


D1  730 


D2  730 


D4  730 


Time(sec) 


Figure  4.  3-channel  NIRS  results  measured  on  a  tumor  dynamic  phantom  with  two- 
size  tubing  networks.  Three  traces  represent  the  readings  at  D1  (in  transmission 
mode  and  located  between  the  two  tube  networks),  D2  (near  the  large  tube),  and  D4 
(near  the  small  tube). 


Figure  4  shows  the  experimental  results  from  another  phantom,  where  the  two  tubes  have 
the  same  diameter  and  length,  but  with  different  controlled  flow  rates.  By  alternating  the  flow 
rate,  we  observed  the  results  similar  to  those  shown  in  Figure  3.  The  first  step  for  this 
experiment  was  to  inject  the  ink  solution  only  into  the  top  tube  (near  D4,  Fig.2)  with  a  flow  rate 
of  20  ml/hr.  As  a  result,  D4  (green  trace)  showed  the  largest  increase  of  Apa,  while  D2  (pink 
trace)  showed  the  smallest  increase  of  Apa  since  D2  is  quite  away  from  the  top  network.  For  the 
second  step,  the  bottom  network  (near  D2)  was  injected  with  the  ink  solution  at  the  same  flow 
rate  of  20  ml/hr.  Here,  D2  detected  the  largest  increase  of  Apa,  and  D4  showed  the  smallest 
increase  of  Apa.  In  both  of  the  above  cases,  the  signals  from  D1  (blue  trace)  showed  similar 
profiles  as  the  others,  without  clear  bi-exponential  patterns.  Then,  as  the  third  step,  the  ink 
solution  was  injected  into  both  of  the  tube  networks  with  the  same  flow  rate  (20  ml/hr). 
Similarly,  we  do  not  observe  any  clear  features  of  the  two-exponential  increase  of  Apa.  For  the 
fourth  step,  we  injected  the  solution  into  the  tubes  with  two  different  flow  rates:  20  ml/hr  for  the 
top  tube  and  5  ml/hr  for  the  bottom  tube.  Now  the  data  clearly  show  a  bi-exponential  behavior 
detected  by  Dl. 


Further  time  constant  analysis  helps  understand  the  meaning  of  the  two-exponential 
feature  measured  at  Dl.  Fig.  5a  shows  a  one-exponential  fitting  of  the  fast  increase  of  Apa 
obtained  at  D4,  with  a  time  constant,  x,  of  20.9  ±  0.5  sec.  Fig.  5b  is  a  similar  fitting  for  the  signal 
at  D2,  with  T=1 2 1  ±3  sec.  Finally,  Fig.  5c  shows  a  two-exponential  fitting  for  the  signal  at  Dl 
with  xl=18.1+  0.9  sec  and  x2=134  ±11  sec  for  both  the  fast  and  slow  components,  respectively. 
The  data  clearly  demonstrate  that  the  time  constant  for  the  fast  component,  xl,  from  Dl  (18,1 
sec)  matches  well  with  x  from  D4  (20.9  sec),  and  x2  from  Dl  (134  sec)  matches  well  with  x  from 
D2  (121  sec).  This  suggests  that  the  dynamic  signals  obtained  from  Dl  with  a  fast  and  slow 
component  indeed  result  from  the  dynamic  changes  in  both  of  the  tubes. 


Detector  #4  (09/27/02)  Detector  #2  (09/27/02) 


|  y=  m1*(1-exp(-m0/m2))  j 

Value 

Error 

ml 

0,25587 

0.0021478 

m2 

120.56 

2.7574 

Chisq 

0,044656 

NA 

R 

0,98517 

NA 

j  y  =  m1*(1-exp(-m0/m2))  ] 

Value 

Error 

ml 

0,17366 

0.00053177 

m2 

20,919 

0.52579 

Chisq 

0,019571 

NA 

R 

0,96538 

NA 

Figure  5  Temporal  profiles  of  the  NIRS  measurements  from  D4  (a)  and  D2  (b)  with 
the  one-exponential  fitting  parameters  shown  at  the  bottom  of  each  figure. 
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Detector  #1  (09/27/02) 


|  y  =  m1*{1-exp(-m0/m2)}+m3*{1,„  | 

Value 

Error 

ml 

0.068109 

0,0024093 

m2 

18.129 

0.87637 

m3 

0.046865 

0.0017792 

m4 

133,65 

11354 

Chisq 

0.0018442 

NA 

R 

0.99287 

NA 

Figure  5c  Temporal  profile  of  the  NIRS  measurement  from  D1 
with  the  bi-exponential  fitting  parameters  shown  on  the  left. 


As  shown  above,  these  dynamic  fast-  and  slow-changes  simulate  very  well  the  changes 
observed  in  our  breast  tumor  oxygenation  measurements,  supporting  our  hypothesis1  that  the  bi- 
phasic  tumor  hemodynamic  feature  results  from  a  well-perfused  and  poorly  perfused  region  in 
the  tumor  vasculature  under  carbogen  interventions.  In  this  way,  the  multi-channel  NIRS  system 
allows  us  to  study  the  dynamic  heterogeneity  at  different  locations,  showing  that  the  tumor 
vasculature  in  the  poorly  perfused  region  (with  a  lower  perfusion  rate)  can  be  eventually 
oxygenated  if  the  two  dynamic  components  are  observed.  On  the  other  hand,  if  only  the  fast 
component  exists,  it  suggests  that  the  poorly  perfused  region  in  the  tumor  may  not  have  much 
chance  to  be  improved  for  its  vascular  oxygenation  under  carbogen  intervention. 


2.  Computer  Simulations  to  Understand  the  Phantom  Experiments 

To  further  support  our  tumor  hemodynamic  phantom  study,  we  used  FEMLAB  software 
to  simulate  the  dynamic  process  of  light  distribution  as  well  as  the  NIR  measurement  in  tissues, 
such  as  in  breast  tumors.  In  this  way,  we  can  better  understand  the  effects  of  tumor  heterogeneity 
on  the  multi-channel  NIR  measurements. 

2.1  Introduction  to  FEMLAB 


Based  on  the  Finite  Element  Method  (FEM),  FEMLAB  is  interactive  computer  software 
and  developed  by  the  COMSOL  Group  (COMSOL  Inc.  Burlington,  MA).  It  uses  numerical 
approaches  to  solve  Partial  Differential  Equations  (PDE)  in  modeling  and  simulating  various 
engineering  problems.  The  FEMLAB  Graphical  User  Interface  (GUI)  is  used  for  the  modeling 
process  by  using  ready-to-use  application  modules  or  by  using  equation-based  modeling.  In 
general,  several  of  ready-to-use  application  modules  are  commercially  available  for  engineering 
problems:  such  as  chemical  engineering  modules,  electromagnetic  module,  and  structural 
mechanics  module.  FEMLAB  solves  different  systems  of  equations  simultaneously  by  using 
multi-physics  modeling,  and  FEMLAB  models  can  also  be  used  with  various  MATLAB 
toolboxes. 
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2.2  Simulation  of  Light  Propagation  in  Tumor  Tissue 

According  to  diffusion  theory,  the  simulation  problem  starts  with  a  PDE  to  describe  the 
light  propagation  in  the  tumor  tissue  as  follows: 

V2W(r)~V(r)  =  S0(r) 

where  \jf(r)  is  the  light  distribution  at  position  r  inside  the  tumor  tissue,  jAa  =  0.1  cm'1  is  the 
absorption  coefficient,  D  =  0.033  cm  is  the  diffusion  coefficient,  and  5'0(r)is  the  light  source 
term.  The  boundary  condition  for  light  propagation  in  the  simulated  tumor  tissue  is  given  by 

y/(r)  -  2  AD  n  .V  =  0 . 

In  the  simulation,  the  Neumann  Boundary  condition  was  used  with  the  various  parameters 
specified  as:  Internal  reflection  A  =  2.9193  and  D  =  0.033  cm.  The  simulation  model  was  then 
meshed  in  the  mesh  mode  using  1 147  elements  and  609  nodes.  Finally,  it  is  solved  using  the 
stationary  nonlinear  solver  type. 

Specifically,  the  geometry  for  the  FEM  tumor  simulations  is  shown  in  Figure  6,  and  it 
mimics  a  2-dimensional  horizontal  cross  section  of  our  tumor  phantom.  El  represents  the  body 
of  tumor  phantom  (diameter  =  4  cm),  and  R5  shows  location  of  the  light  source.  Furthermore, 
each  tube  used  in  the  dynamic  phantom  is  represented  by  two  rectangles  (0.1  x  2,4  cm  rectangle) 
as  a  blood  vessel  network:  R1  and  R2  are  representing  vessels  with  a  fast  flow  rate,  and  R3  and 
R4  denote  vessels  with  a  slow  flow  rate. 


Figure  6  Geometry  for  FEM  simulations  of  a  breast  tumor. 

The  diffusion  equation  has  been  selected  to  predict  photon  distribution  in  the  phantom  as 
the  ink  solution  flows  through  the  tubes.  Absorption  coefficients  for  the  background  phantom 
and  tubes  were  0.03  and  1.5  cm'1,  respectively.  We  also  used  0.03  cm  as  the  diffusion  coefficient 
of  the  phantom  since  it  gives  a  reduced  scattering  coefficient  of  10  cm'1. 
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2.3  Simulation  Results  of  the  Tumor  Dynamics 


Different  flow  rates  used  in  the  phantom  study  were  simulated  by  changing  the  tube  size 
with  different  frame  rates.  In  other  words,  a  fast  flow  rate  was  simulated  by  increasing  the  tube 
length  with  a  large  increment  per  frame  (e.g.,  0  to  2.4  cm  with  0.6  cm  increment  per  frame),  and 
a  slow  flow  rate  was  simulated  by  increasing  tube  length  with  a  small  increment  (0  to  2.4  cm 
with  0.02  cm  increment  per  frame).  In  this  way,  a  flow  rate  difference  between  the  two  tubes  can 
be  established. 

Figure  7  shows  an  example  of  the  FEM  tumor  phantom  simulation.  Figure  7a  is  the  light 
distribution  of  a  dynamic  phantom,  where  there  was  ink  injection  into  only  the  left  tube  with  a 
fast  flow  rate.  Figure  7b  is  the  simulation  when  the  ink  solution  flowed  into  only  the  right  tube 
with  a  slow  flow  rate.  With  the  same  frame  rate,  Figure  7a  is  the  result  of  the  simulation  at  the 
6th  frame,  while  Figure  7b  is  simulated  for  the  120th  frame.  This  implies  that  the  left  tube  has  a 
flow  rate  20  times  faster  than  the  right  tube. 


Figure  7  FEM  simulations  of  light  distribution  inside  of  a  dynamic  phantom.(a):  the 
results  with  a  fast  simulated  flow  rate;  (b):  the  results  with  a  slow  flow  rate  at  a  later  time. 


Figure  8  FEM  simulations  of  light  distribution  inside  of  a  dynamic  phantom. 

(a):  the  results  with  a  fast  simulated  flow  rate;  (b):  the  results  with  a  slow  flow  rate. 
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We  also  changed  the  position  of  the  light  source  to  see  the  effect  of  vessel/tube  geometry 
on  bi-phasic  feature  of  the  tumor  hemodynamic  measurements.  We  changed  the  position  of  the 
light  source  to  the  right  side  of  the  phantom  (Fig.  8a)  and  to  the  center  of  the  phantom  (Fig.  8b). 
In  the  geometry  of  Fig.  8a,  the  light  penetrates  the  slow-flow  object  first  and  then  the  fast-flow 
object.  This  simulation  was  done  to  find  whether  in  this  geometry,  the  light  signal  changes  still 
have  the  two-exponential  behavior.  In  other  words,  we  wish  to  investigate  whether  the  bi-phasic 
hemodynamic  feature  depends  on  the  orientations  of  the  vessels/tube  networks.  We  are  in  the 
continuation  for  this  study.  Furthermore,  we  are  also  in  the  process  of  comparing  the  simulation 
results  with  those  obtained  in  the  phantom  experiments,  particularly  to  compare  the  time 
constants. 


3.  Simultaneous  Measurement  using  NIR  and  19F  EPI  relaxation  of  HFB 

To  prepare  ourselves  for  Aim  4,  which  is  to  study  the  influence  of  three  interventions  on 
Hb02  and  p02  of  the  tumor  simultaneously,  we  have  conducted  the  following  work  for 
simultaneous  measurements  using  multi-channel  NIRS  and  19F  EPI  relaxation  of  HFB 
(Hexafluorobenzene). 

3.1  Animal  &  tumor  model: 


NF13762  mammary  adenocarcinomas  were  implanted  in  skin  pedicles  on  the  forebacks 
of  the  female  Fisher  344  rats  (~150g)  until  tumors  reach  approximately  1  cm  in  diameter.  The 
measurements  would  be  taken  when  the  tumor  had  a  diameter  of  1  cm,  1.5  cm,  and  2  cm.  Rats 
were  anesthetized  with  150  jxL  Ketamine  hydrochloride  (100  mg/ml,  i.p.)  and  maintained  under 
general  gaseous  anesthesia  with  1.3%  isoflurane  in  air  (1  dm3/min). 

20  jiL  of  Hexafluorobenzene  (PCR  Inc.,  Fainesville,  FL)  was  injected  directly  into  the 
tumor,  either  centrally  or  peripherally  in  a  single  plane  using  a  Hamilton  syringe  (Reno,  NV) 
with  a  sharp  needle  (32  gauge).  The  needle  was  inserted  manually  to  penetrate  across  the  whole 
tumor  and  withdrawn  ~1  mm  to  reduce  the  tissue  pressure,  and  2~5  pL  of  HFB  were  injected. 
The  needle  was  repeatedly  withdrawn  at  a  step  size  of  1~2  mm  and  further  HFB  injected.  A  total 
of  2~3  tracks  of  HFB  were  injected  in  the  form  of  a  fan  shape. 

The  rat  was  placed  on  its  side  in  a  cradle  with  a  thermal  blanket  to  maintain  body 
temperature  at  37°C.  A  home-built  tunable  and  matchable  2  cm  cooper  coil  was  made  and  placed 
around  the  tumor  with  the  tumor  centered  inside  the  coil.  In  addition,  a  sealed  capillary  tube  of 
HFB  doped  with  benzene  (as  a  shift  agent)  was  placed  alongside  the  tumor.  This  provided  a 
quantitative  standard  with  approximately  3  ppm  offset  from  the  intra-tumoral  HFB,  allowing 
clearance  from  the  tumor  to  be  assessed.  A  home-built  plastic  NIR  source-detector  fixer  was  also 
made  to  fix  the  two  NIR  fiber  bundle  tips  around  the  receiver  coil  within  the  MRI  facility,  as 
shown  in  Figure  9. 
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Fig  9,  Home-built  plastic  NIR  source-detector  fixer  with  a  receiver  coil, 

3.2  MRI  scanner  and  RF  coils 


All  in  vivo  MRI  experiments  were  performed  on  an  Omega  CSI  4,7  T  superconducting 
magnet  with  actively  shielded  gradients.  The  magnet  operated  at  a  proton  (*H)  frequency  of 
200.106  MHz  and  a  fluorine  (19F)  frequency  of  188.273  MHz.  It  consisted  of  the  Omega  software 
that  can  accomplish  data  acquisition  and  analysis  with  the  UNIX  operating  system.  Two  single- 
turn  solenoid  coils  (1.5  cm  and  2,5  cm)  were  built  for  this  study.  Their  sizes  matched  the  tumor 
sizes. 


3.3  NIR  spectroscopy 

NIR  spectroscopy  uses  commercially  available  in-phase  and  quadrature  (IQ)  demodulator 
chips  to  demodulate  the  detected,  amplitude-modulated  optical  signals.  22-feet  long  fibers  with  a 
diamagnetic  material  for  the  probe  tips  were  used  within  the  magnetic  field  to  transmit  the  laser 
light  to  the  tumor  and  to  detect  the  optical  signal  from  the  tumor.  The  reason  to  use  22 -feet  long 
fibers  is  to  keep  the  NIR  oximeter  away  from  the  magnetic  field.  Three  of  home-built  converters 
have  been  made  for  the  long  fiber  bundles  to  couple  with  the  home-built  connector. 

3.4  Investigation  of  Tumor  pQ?  with  FREDOM  and  HbO?  with  NIR 

Currently  for  this  study,  we  are  in  the  process  of  using  the  FREDOM  approach 
(Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)  based  on 
hexafluorobenzene  (HFB)  to  measure  tumor  p02  and  the  NIRS  system  to  measure  Hb02 
simultaneously. 

The  respiratory  challenge  paradigm  for  the  rats  will  be: 

air  ->  100%  O2  air  ->  95%  02  +5%  C02  (carbogen)  ->  air 
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Key  Research  Accomplishments 

1)  We  have  obtained  and  tested  a  multi-channel  NIRS  system  for  simultaneous  dynamic 
oxygenation  measurements  at  different  locations  of  breast  tumors. 

2)  We  have  developed  breast  tumor  hemodynamic  phantoms  and  conducted  multi-channel 
NIRS  measurements  on  the  newly  developed  phantoms  to  simulate  the  dynamic  HbC>2 
changes  of  rat  breast  tumors  under  hyperoxic  inhalations.  This  study  supports  our  earlier 
dynamic  model  to  explain  the  bi-phasic  dynamic  oxygenation  process  of  breast  tumors 
under  hyperoxic  interventions. 

3)  We  have  further  implemented  and  utilized  FEMLAB  computer  software  to  mimic 
dynamic  light  distribution  in  a  simulated  tissue  vasculature  so  as  to  understand  the  effects 
of  tumor  heterogeneity  on  the  multi-channel  NIR  measurements. 

4)  We  have  implemented  home-built  tunable  2  cm  cooper  coil,  which  will  be  placed  around 
the  rat  breast  tumor  within  the  MRI  facility  during  the  simultaneous  measurements 
between  FREDOM  and  NIRS.  Also,  a  home-built  plastic  NIR  source-detector  coupler 
(fixer)  was  made  to  hold  the  two  NIR  fiber  bundle  tips  around  the  receiver  coil  within  the 
MRI.  This  work  is  done  to  prepare  for  experiments  needed  for  Aim  4. 
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Conclusions 

From  the  work  that  we  have  conducted  in  Year  2,  we  can  draw  the  following  conclusions: 

1)  With  a  multi-channel  NIRS  system  and  the  newly  developed  tumor  hemodynamic 
phantoms,  we  have  experimentally  demonstrated  that  the  previously  observed,  bi-phasic 
temporal  changes  in  Hb02  of  breast  tumors  under  carbogen  intervention  can  be  associated 
with  two  perfusion  rates  or  two  blood  flow  rates. 

2)  The  newly  developed  hemodynamic  phantoms  are  helpful  to  simulate  the  hemodynamic 
processes  undertaken  within  the  tumor  tissues,  and  the  phantoms  will  be  further  used  in  Year 
3  study. 

3)  FEMLAB  has  proven  to  be  valuable  in  simulating  the  light  propagation  in  the  breast 
tumor  tissue  and  the  corresponding  NIRS  measurements,  and  we  plan  to  use  FEMLAB  for 
further  experimental  validations  in  both  laboratory  and  animal  measurements. 

4)  After  the  implementation  of  the  home-built,  tunable  2-cm  cooper  coil  and  a  home-built 
plastic  NIR  source-detector  coupler,  it  is  feasible  and  ready  to  conduct  the  simultaneous 
measurements  between  FREDOM  and  NIRS. 
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Abstract,  This  study  investigates  the  correlation  of  tumor  blood  oxy¬ 
genation  and  tumor  p02  with  respect  to  carbogen  inhalation.  After 
having  refined  and  validated  the  algorithms  for  calculating  hemoglo¬ 
bin  concentrations,  we  used  near-infrared  spectroscopy  (NIRS)  to 
measure  changes  of  oxygenated  hemoglobin  concentration 
(A[Hb02])  and  used  an  oxygen  needle  electrode  and  19F  MRI  for 
p02  measurements  in  tumors.  The  measurements  were  taken  from 
Dunning  prostate  R3327  tumors  implanted  in  rats,  while  the  anesthe¬ 
tized  rats  breathed  air  or  carbogen.  The  NIRS  results  from  tumor  mea¬ 
surements  showed  significant  changes  in  tumor  vascular  oxygenation 
in  response  to  carbogen  inhalation,  while  the  p02  electrode  results 
showed  an  apparent  heterogeneity  for  tumor  pC)2  response  to  carbo¬ 
gen  inhalation,  which  was  also  confirmed  by  19F  MR  p02  mapping. 
Furthermore,  we  developed  algorithms  to  estimate  hemoglobin  oxy¬ 
gen  saturation,  s02,  during  gas  intervention  based  on  the  measured 
values  of  A[Hb02]  and  p02.  The  algorithms  have  been  validated 
through  a  tissue-simulating  phantom  and  used  to  estimate  the  values 
of  s02  in  the  animal  tumor  measurement  based  on  the  NIRS  and 
global  mean  p02  values.  This  study  demonstrates  that  the  NIRS  tech¬ 
nology  can  provide  an  efficient,  real-time,  noninvasive  approach  to 
monitoring  tumor  physiology  and  is  complementary  to  other  tech¬ 
niques,  while  it  also  demonstrates  the  need  for  an  NIR  imaging  tech¬ 
nique  to  study  spatial  heterogeneity  of  tumor  vasculature  under  thera¬ 
peutic  interventions,  ©  2003  Society  of  Photo-Optical  Instrumentation  Engineers . 
(DOl:  10.1117/1.15270491 
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1  Introduction 

It  has  long  been  known  that  hypoxic  tumor  cells  are  more 
resistant  to  radiation  therapy  than  well -oxygenated  tumor 
cells.1  Breathing  elevated  oxygen  (100%)  or  carbogen  (95% 
02 ,  5%  C02)  has  been  used  during  therapy  for  an  attempt  to 
improve  tumor  oxygenation.2*3  To  monitor  tumor  tissue  oxy¬ 
gen  tension4  and  its  dynamic  changes  under  respiratory  inter¬ 
ventions,  various  methods  are  available,  including  fiber  optic 
sensors,5  oxygen  electrodes,6  and  electron  spin  resonance.7 
MRI  has  the  further  advantage  of  providing  dynamic  maps  of 
p02,  which  can  reveal  tumor  heterogeneity.8  While  NIRS 
does  not  quantify  p02,  it  can  indicate  dynamic  changes  in 
vascular  oxygenation  and  has  the  advantage  of  being  entirely 
noninvasive,  providing  real-time  measurements,  and  being 
cost-effective  and  portable.  Furthermore,  it  would  be  impor¬ 
tant  to  correlate  the  changes  between  tissue  p02  and  vascular 
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oxygenation  of  the  tumors  since  little  is  known  about  oxygen 
transfer  from  the  tumor  vasculature  to  tumor  tissue. 

The  basic  principle  of  NIRS  rests  on  the  fact  that  oxygen¬ 
ated  and  deoxygenated  hemoglobin  molecules  are  major  chro- 
mophores  in  tissue  in  the  near-infrared  region  (700  to  900 
nm),  and  they  exhibit  distinct  absorption  characteristics.  In 
principle,  the  concentrations  of  oxygenated  hemoglobin 
[Hb02],  deoxygenated  hemoglobin  [Hb],  and  oxygen  satura¬ 
tion  of  hemoglobin  s02  can  be  determined  by  measuring  light 
absorption  and  scattering  in  tissue  based  on  diffiision  theory. 
However,  the  theory  works  well  only  for  large  and  homoge¬ 
neous  media,9fl°  Accurate  quantification  of  tumor  oxygenation 
in  our  approach  is  currently  limited  to  relative  changes  in 
[Hb02]  and  [Hb]  due  to  considerable  heterogeneity  and  finite 
size  of  tumors. 

The  goal  of  this  study  was  to  investigate  the  correlation  of 
tumor  blood  oxygenation  and  tumor  p02  in  response  to  car- 
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bogen  intervention  and  to  develop  a  suitable  algorithm  to  es¬ 
timate  the  hemoglobin  oxygen  saturation  of  the  tumor  under 
the  intervention.  Specifically,  in  Sec,  2  of  this  paper,  we  de¬ 
rive  accurate  expressions  for  calculating  changes  in  [Hb02] 
and  [Hb]  to  compensate  for  the  differences  in  optical  path 
length  at  two  wavelengths  and  an  algorithm  to  estimate  abso¬ 
lute  sO?  values  of  the  tumor  during  gas  intervention.  The 
algorithms  are  validated  through  tissue-simulating  phantoms 
and  used  to  estimate  tumor  s02  in  the  animal  measurement 
using  the  NTRS  and  mean  p02  values,  as  mentioned  in  Secs.  3 
and  4.  In  Sec,  4,  we  will  show  that  while  NIRS  results  tended 
to  be  similar  for  several  tumors,  p02  electrode  measurements 
showed  considerable  variation  even  in  the  same  tumor  type, 
suggesting  distinct  tumor  heterogeneity.  In  Sec,  5,  we  discuss 
the  need  to  develop  an  NIR  imaging  technique  in  order  to 
study  spatial  heterogeneity  of  tumor  vasculature  under  oxygen 
interventions.  Finally,  we  conclude  that  the  NIRS  technology 
can  provide  an  efficient,  real-time,  noninvasive  approach  to 
monitoring  tumor  physiology  and  is  complementary  to  other 
techniques. 


2  Theory  and  Algorithm  Development 

2,1  Algorithms  to  Quantify  Changes  in  [Hb02]  and 
[Hb] 

NIR  spectroscopy  can  be  used  to  measure  hemoglobin  con¬ 
centrations  and  oxygen  saturation  since  light  absorptions  of 
Hb02  and  Hb  are  different  at  the  wavelengths  selected  (75B 
and  785  nm).  In  common  with  our  previous  work,11  we  as¬ 
sumed  that  Hb02  and  Hb  are  the  only  significant  absorbing 
materials  in  tumors  within  the  selected  NIR  range  of  700  to 
900  nm.  Based  on  Becr-Lambert’s  law,  the  absorption  coeffi¬ 
cients  jxa  comprise  the  extinction  coefficients  for  deoxyhemo¬ 
globin  (eHb)  and  oxyhemoglobin  (£Hb02.)  multiplied  by  their 
respective  concentrations: 

^5S=  2.3{e™[Hb]  +  45b802[Hb02]},  (1) 

nT=  2.3{S?[Hb]  +  e7t502[Hb02]},  (2) 

where  the  factor  of  2.3  results  from  the  different  definitions  of 
fAa  and  e  in  relation  to  the  incident  and  detected  optical  in¬ 
tensities.  The  conventional  definitions  for  fxa  and  e  are  I 
—  I0  exp(—  ix JS)  and/=/o10“ECi%  respectively,  where  /0  and 
I  are  the  incident  and  detected  optical  intensities  in  transmis¬ 
sion  measurement  of  a  nonscattering  medium,  C  is  the  con¬ 
centration  of  hemoglobin  measured  in  millimoles  per  liter, 
and  L  is  the  optical  path  length  through  the  medium  in  centi¬ 
meters.  Therefore,  we  should  have  a  relationship  of  f ia 
=  2.3  eC. 

We  have  not  yet  completed  a  suitable  algorithm  to  com¬ 
pute  of  rat  tumors  due  to  their  finite  size  and  high  hetero¬ 
geneity,  Instead  of  diffusion  theory,  we  modified  Beer- 
Lamberfs  law,  i.e.,  /ifl  =  2,3eC=(2,3/L)!og(/0//),  to  analyze 
the  data  using  only  amplitude  values  to  quantify  changes  in 
[Hb02]  and  [Hb],  In  this  case,  /0  is  the  detected  light  inten¬ 
sity  when  no  absorption  is  present  Specifically,  changes  in 
absorption  coefficient  of  the  tumor,  A  /j,a ,  between  baseline 
and  transient  conditions  under  respiratory  intervention  can  be 
expressed  as 


AM0=/i0r_M0fl=2.31og(^5//(r)/Z,,  (3) 

where  L  is  the  optical  path  length  and  A  B  and  A  T  are  baseline 
and  transient  amplitudes  of  the  measured  optical  signals,  re¬ 
spectively. 

By  manipulating  Eqs.  (1)  to  (3),  changes  of  [Hb02]  and 
[Hb]  due  to  an  intervention  can  be  expressed  using  the  trans¬ 
mitted  amplitudes  of  the  light  through  the  tumor  as: 


A[Hb02]  =  —  11,73* 


log  (Ab/At)™ 

£758 


+ 14.97* 


log  (AB/AT)m 


r  785 


(4) 


r  ,  log(  AgIAT)m  log(  Ab/At)7SS 

A[Hb]  =  8.09*  75S - 6.73*- 


r  7S5 


(5) 


where  Z,758  and  I785  are  optical  path  lengths  between  the 
source  and  detector  at  758  and  785  nm,  respectively.  The  units 
of  A[Hb02]  and  A[Hb]  in  Eqs.  (4)  and  (5)  are  in  millimolar. 
The  constants  given  in  the  equations  were  computed  with  the 
extinction  coefficients  for  oxygenated  and  deoxygenated  he¬ 
moglobin  at  the  two  wavelengths  used.12  The  constant  values 
are  slightly  different  from  our  previous  report11  due  to  a  slight 
shift  in  wavelength  (782  to  785  nm)  from  one  laser  source, 
but  the  actual  differences  between  the  values  of  A[HbOa]  and 
A[Hb]  calculated  from  our  previous  report  and  from  Eqs,  (4) 
and  (5)  are  little  and  negligible. 

In  principle,  L758  and  L785  given  in  Eqs,  (4)  and  (5)  are  not 
constants,  depending  on  both  the  source-detector  separation 
and  the  optical  properties  of  the  measured  medium.  Optical 
path  length  in  a  scattering  medium  L  has  been  expressed13  as 
the  source-detector  separation  d  multiplied  by  a  differential 
pathlength  factor  (DPF),  i.e.,  L  =  d*DPF ,  DPF  values  of 
blood-perfused  tissues  should  be  wavelength-  and 
oxygenation-dependent,  and  they  have  been  studied  inten¬ 
sively  for  muscles14  and  brains15  with  approximate  values  of  4 
to  6  and  5  to  6,  respectively  Little  is  known  about  DPF  for 
tumors  although  a  DPF  value  of  2,5  has  been  used  by  others.16 
In  our  approach,  we  define  two  parameters,  /3Hb02  Pm  ? 

as  ratios  between  DPF1S%  and  DPF7BS  for  oxygenated  blood 
and  deoxygenated  blood,  respectively,  as  given  below: 


DPF 758 
Arb02-I  '^ppm 


Hb02 


I758 

JIB 


Hb02 


Pm~ 


758 


DPF 

~5¥Fm 


Hb 


L758 

jim 


Hb 


(6) 


Substituting  Eq.  (6)  into  Eqs,  (4)  and  (5)  leads  to 


A[Hb02]=- 


11.73 

PubQ2 


758 


logf  —  J  +14.97  log 


\  785 


Ar) 


d X DPF a 


(7) 
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A[Hb]= 


8.09  (Ab\15S 

^'0Sfe/ 


6.73  log^  — 


785 


dXDPFft 


(8) 


where  DPFQ  is  a  mean  DPF  at  785  nm  for  both  oxygenated 
and  deoxygenated  states,  i.e.,  DPF0=DPF7^O2=DPF7^ , 
which  is  assumed  to  be  the  same  for  both  A[Hb02]  and 
A[Hb].  This  assumption  is  based  on  the  fact  that  the  absorp¬ 
tion  difference  between  oxygenated  and  deoxygenated  blood 
at  785  nm  is  much  smaller  than  that  at  758  nm.  The  maximal 
relative  error  caused  by  this  assumption  in  tumor  oxygen  in¬ 
terventions  was  estimated  to  be  less  than  12%,  and  detailed 
justification  and  discussion  were  given  in  Ref  11.  Since  our 
focus  is  on  dynamic  changes  in  tumor  [Hb02]  under  carbo- 
gen  intervention,  we  simplify  Eqs.  (7)  and  (8)  to  Eqs,  (9)  and 
(10)  by  including  DPF  §  in  the  unit: 


11.73  M/758 


A[Hb02]=  ■ 


Pmoi  \  A  t 


log  -r  + 14.97  lof 
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(10) 


A[HbtotaI]  =  A[Hb02]  +  A[Hb] 

-uM^r+8-24i°sfer 

d 

(15) 

Equations  (13)  to  (15)  will  be  used  in  calculating  A[Hb02], 
A[Hb],  and  A[Hbtota|]  in  tissue  phantoms  and  tumors  during 
gas  interventions  in  this  paper. 

The  units  for  A[Hb02],  A[Hb],  and  A[Hbtotal]  in  Eqs. 
(13)  to  (15)  are  mM/DFF0,  which  is  still  a  variable,  depend¬ 
ing  on  the  optical  properties  of  the  tumor  at  a  particular  wave¬ 
length.  Since  our  study  involves  changes  in  [Hb02]  due  to 
respiratory  challenges,  we  can  obtain  a  normalized  A[Hb02] 
at  its  maximal  value,  i.e.,  A[Hb02]/A[Hb02]max,  to  elimi¬ 
nate  the  unit  so  as  to  minimize  the  effect  of  DPF  on  our 
results.  Next,  we  will  show  that  a  normalized  A[Hb02]  has  a 
close  relationship  with  hemoglobin  oxygen  saturation  s02 . 


2*2  Relationship  Among  Normalized  A[Hb02],  s02 
and  Blood  p02 

We  define  s02  values  of  the  measured  sample  at  the  baseline, 
transient  state,  and  maximal  state,  i.e.,  (s02)base,  (s02)/ ,  and 
(s02)max ,  respectively: 


where  the  units  for  Eqs.  (9)  and  (10)  become  mM/DPF0. 

To  further  quantify  /?Hb02  Pm  >  we  associate  L  to  pba 
by  L  =  (s/3f2)d  {fifJ ^)I/2,  where  ft'  is  the  reduced  scatter¬ 
ing  coefficient,  according  to  Sevick  et  al.10  and  Liu.!?  Equa¬ 
tion  (6)  becomes 
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where  [xa=23  sC  and  values  at  two  wavelengths  are 
canceled,  assuming  that  /*'  (758  nm)s^  (785  nm).  By 
calculating  the  hemoglobin  extinction  coefficients  at  758  and 
785  nm,12  we  obtained  ^^02=  1.103  and  ^=0.9035.  Sub¬ 
stituting  these  values  into  Eqs.  (9)  and  (10)  results  in  the  final 
expressions  for  A[Hb02]  and  A[Hb]: 


A[Hb02]  =  - 


AR\m  IAb 


10.63  log  -H  +  14.97  log \-f- 
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(13) 


A[Hb]=  ■ 
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8.95  log]  —  j  —  6.73  log|  —  j 
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■■  04) 


A[Hbtoai]  can  also  be  obtained  by  adding  Eqs.  (13)  and  (14): 


(sC^max- 


[HbQ2]max 

[Hbtotal] 

max 


(18) 


where  [Hb02]bases  (jHb02]r  5  and  [Hb02]max  are  correspond¬ 
ing  to  oxygenated  hemoglobin  concentrations  at  the  respec¬ 
tive  state.  Mathematically,  it  follows  that 


As°2  ^  (s02),-(s02)basc 

^  s^2  max  (  ®^2  )  max  -  (s02)base 

/  [HbQ2],  [HbQ2]base\ 

__  \  [blbtotai]r  [Hbtotallb  as  e/  .... 

~  /  [HbQ2]max  [HbQ2]  base  \ 

\  [Hbtotal]  max  [Hbtota1]base/ 

During  a  cycle  of  oxygenation  and  deoxygenation  in  a 
blood-perfused  tissue,  if  the  total  concentration  of  hemoglobin 
remains  constant,  we  have  the  following  condition: 
[ Hbtota1] max -  [Hbtota  J t  -  [ Hbtotal] basc .  In  the  case  of  tumors 
under  gas  intervention,  total  hemoglobin  concentration  does 
not  always  remain  constant,  but  the  changes  in  [Hb]totai  ap¬ 
peared  relatively  small  in  comparison  to  the  changes  in 
[Hb02]. 11,18  It  is  reasonable  to  assume  that  A[HbtotaJ 
^[HbtotaI],  i.e.,  the  condition  of  [Hbtotal]max=fHbtotai]f 
=  [HbtotaJ]basc  still  holds  for  the  tumor  under  oxygen/carbogen 
interventions.  Then,  Eq.  (19)  becomes 
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AsQ2  (s02) -(s02)tee  A[HbQ2] 

^  SO2  max  (  sOj  )  max  ~  (  s^2  )  base  A  [  Hb02]  max 

(20) 

To  further  make  correlation  between  the  normalized 
A[Hb02]  and  blood  p02,  HilTs  equation19  can  be  combined 
with  Eq.  (20)  to  characterize  oxygen  transport  in  the  tissue 
vasculature: 

(pQfr  .  n . 

A[HbQ2]  (Pg0)B+(pOf)”  2base 
A[Hb02]max  —  (s02)max-(s02)base 

(ppfr 

_  (Pf0)"+(pOf)» 

a-b  ’  (Zl) 

where  pOf  is  the  oxygen  partial  pressure  in  blood,  Pf0  is  the 
oxygen  partial  pressure  in  blood  at  sO2  =  50%,  n  is  the  Hill 
coefficient,  a  —  (s02)jnax ,  and  b  —  (s02)basc ,  This  equation  as¬ 
sociates  the  normalized  AHb02  to  blood  p02  in  tissues.  It 
indicates  that  normalized  A[Hb02]  measured  from  tissues/ 
tumors  under  gas  interventions  is  associated  with  normalized 
$02  between  (s02)base  and  (s02)max  of  the  tissue/tumor,  and  it 
predicts  the  relationship  between  the  normalized  A[Hb02] 
and  blood  p02  values  in  the  tissue/tumor  vasculature. 

In  our  phantom  studies,  the  measured  p02  values  are  con¬ 
sidered  as  blood  p02  in  tissue  vasculature  since  blood  is  well 
mixed  in  the  solution  (see  details  in  Sec.  3.4).  Therefore,  val¬ 
ues  of  Pf0,  n}  a ,  and  b  in  Eq.  (21)  can  be  fitted  to  the  experi¬ 
mental  data,  allowing  us  to  determine  the  initial,  transient,  and 
maximal  values  of  s02  of  the  simulating  tissue  due  to  oxygen/ 
nitrogen  interventions, 

2,3  Relationship  Between  Normalized  A[Hb02]  and 
Tissue/Tumor  p02 

In  principle,  blood  p02  and  tissue  p02  are  different,  depend¬ 
ing  on  the  relative  distance  between  a  capillary  vessel,  oxygen 
consumption,  and  the  location  where  p02  is  measured.19  It  is 
shown  that  there  exists  a  constant  pressure  drop  between 
blood  pO?  and  tissue  p02  as  the  blood  passes  through  a  cap¬ 
illary  vessel  Therefore,  it  is  reasonable  to  assume 

pOf=ar-pOj,  (22) 

where  pOf  and  pO[  are  blood  p02  and  tissue  p02  values, 
respectively,  and  a  is  a  constant  representing  an  oxygen  par¬ 
tial  pressure  drop  from  blood  p02  to  a  local  tissue  p02 .  Sub¬ 
stituting  Eq.  (22)  in  Eq,  (21)  results  in 

(pop" 

A[Hb02]  (P^  +  CpOf)'1 

A[Hb02]max  a-b  ’  (23) 

where  P5Q  is  the  oxygen  partial  pressure  in  tissue  resulting 
from  Pf0 ,  and  the  meanings  of  n,  a,  and  b  remain  the  same  as 
in  Eq.  (21).  This  equation  shows  how  normalized  A[Hb02] 
measured  from  tissue  under  gas  interventions  is  associated 
with  both  tissue  p02  and  normalized  s02  between  (s02)base 
and  (s02)max  in  the  tissue  vasculature. 


Ideally,  when  both  A[Hb02]  and  tissue  p02  are  measured 
at  the  same  physical  location,  the  maximal  and  initial  oxygen 
saturations,  i.e.,  a  and  b  in  Eq.  (23),  of  the  measured  tissue 
vasculature  can  be  obtained  by  fitting  Eq,  (23)  to  the  mea¬ 
sured  data.  In  our  tumor  study,  we  then  can  estimate  the  maxi¬ 
mal  and  initial  hemoglobin  oxygen  saturations  of  the  tumor 
by  fitting  the  measured  values  of  global  normalized  A[Hb02] 
and  global  tissue  p02 ,  which  result  from  adding  up  all  local 
pO?  values  obtained  from  l9F  MR  p02  mapping, 

3  Materials  and  Methods 

3.1  Tumor  Model 

Dunning  prostate  rat  tumors  (eight  R3327-HI  and  four 
R3327-AT1)20  were  implanted  in  pedicles  on  the  forebaek  of 
adult  male  Copenhagen  rats,  as  described  in  detail 
previously,21  Once  the  tumors  reached  approximately  1  cm  in 
diameter,  the  rats  were  anesthetized  with  0,2  ml  ketamine  hy¬ 
drochloride  (100  mg/mL;  Aveco,  Fort  Dodge,  I  A)  and  main¬ 
tained  under  general  gaseous  anesthesia  with  isoflurane  in  air 
(1,3%  isoflurane  at  1  dmVmin  air)  through  a  mask  placed  over 
the  mouth  and  nose.  Tumors  were  shaved  to  improve  optical 
contact  for  transmitting  light.  Body  temperature  was  main¬ 
tained  by  a  warm  water  blanket  and  was  monitored  by  a  rec¬ 
tal  ly  inserted  thermal  probe  connected  to  a  digital  thermom¬ 
eter  (Digi-Sense,  model  91100-50,  Cole-Parmer  Instrument 
Company,  Vernon  Hills,  IL).  A  pulse  oximeter  (model  8600, 
Nonin,  Inc.,  Plymouth,  MN)  was  placed  on  the  hind  foot  to 
monitor  arterial  oxygenation  (Sa02).  Tumor  volume  V  (in 
cubic  centimeters)  was  estimated  as  F=(4tt/3)  [(£  +  IT 
+H)16}3,  where  L,  Wt  and  H  are  the  three  respective  or¬ 
thogonal  dimensions. 

In  general,  the  source-detector  fiber  separation  was  about  1 
to  1.5  cm  In  transmittance  geometry,  and  thus  the  maximal 
tumor  volume  interrogated  by  NIR  light  can  be  estimated  as 
follows.  By  the  diffusion  approximation,  the  optical  penetra¬ 
tion  depth  from  the  central  line  between  the  source  and  detec¬ 
tor  is  about  one  half  of  the  separation  (source-detector 
separation^),  The  total  tumor  volume  interrogated  by  NIR 
light  can  be  estimated  as  the  spherical  volume  with  a  radius  of 
one  half  of  d ,  i.e.,  ( wI6)d 3.  In  this  way,  the  estimated  tumor 
volume  interrogated  by  NIR  light  is  in  the  range  of  0.5  to  2,0 
cm3,  depending  on  the  actual  source-detector  separation. 

3.2  N1RS  and  p02  Needle  Electrode  Measurements 

Figure  1  shows  the  schematic  setup  for  animal  experiments 
using  both  NIRS  and  a  p02  needle  electrode.  A  needle  type 
oxygen  electrode  was  placed  in  the  tumor,  and  the  reference 
electrode  was  placed  rectally.  The  electrodes  were  connected 
to  a  picoammeter  (Chemical  Microsensor,  Diamond  Electro- 
Tech  Inc.,  Ann  Arbor,  MI)  and  polarized  at  -0,75  V,  Linear 
two-point  calibrations  were  performed  with  air  (21%  02)  and 
pure  nitrogen  (0%  02)  saturated  saline  buffer  solutions  before 
the  electrode  was  Inserted  into  the  tumor,  and  we  estimated  an 
instrumental  precision  of  2  to  3  mm  Hg.  Measurement  points 
of  p02  were  manually  recorded,  while  the  NIRS  data  were 
acquired  automatically.  Measurements  of  p02  and  NIRS  were 
initiated,  while  the  rats  breathed  air  for  —10  min  to  demon¬ 
strate  a  stable  baseline.  The  inhaled  gas  was  then  switched  to 
carbogen  for  15  min  and  switched  back  to  air. 
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Fig,  1  Schematic  experimental  setup  of  one-channel,  near-infrared, 
frequency-domain  IQ  instrument  for  tumor  investigation  in  vivo.  The 
5-mm-diameter  fiber  bundles  deliver  the  laser  light,  comprising  two 
wavelengths  (758  and  785  nm),  and  detect  the  laser  light  transmitted 
through  the  implanted  tumor.  The  p02  needle  electrode  measures 
tumor  tissue  p02 . 


Our  NIR  system  as  shown  in  Figure  1  (Refs.  1 1  and  22)  is 
a  homodyne  frequency-domain  photon  migration  system 
(NIM,  Inc.,  Philadelphia,  PA)  and  uses  commercially  avail¬ 
able  in-phase  and  quadrature  (IQ)  demodulator  chips  to  de¬ 
modulate  the  detected,  amplitude-modulated  optical  signal. 

33  Experimental  Validation  for  /?Hb02  and  Pm 
Values 

In  order  to  validate  /?Hb02  and  f3m  values,  we  conducted 
phantom  calibration  measurements.  We  used  2  1  of  0.01  M 
phosphate  buffered  saline  (P-3813,  Sigma,  St  Louis,  MO)  and 
1%  Intralipid  (Intralipid®  20%,  Baxter  Healthcare  Corp., 
Deerfield,  IL)  with  pH—7.4  at  25  °C,  To  deoxygenate  the  so¬ 
lution,  14  g  of  baking  yeast  was  dissolved  in  the  phantom 
solution,  and  pure  oxygen  gas  was  used  to  oxygenate  the  so¬ 
lution.  After  the  yeast  was  well  mixed  in  the  solution,  3  ml  of 
human  blood  was  added  twice.  When  the  blood  was  fully 
deoxygenated,  pure  oxygen  was  introduced  in  the  solution  to 
oxygenate  the  blood.  After  the  blood  was  fully  oxygenated, 
oxygen  blowing  was  stopped  in  order  to  deoxygenate  the  so¬ 
lution  with  yeast  again. 

Equations  (4)  and  (5)  were  applied  to  the  raw  amplitude 
data  to  calculate  A[Hb02]  and  A[Hb].  Large  unexpected  and 
erroneous  fluctuation  of  A[Hbtotal]  (= A[Hb02]  +  A[Hb]) 
were  seen  during  the  oxygenation  and  deoxygenation  cycles 
(Figure  2).  However,  when  we  applied  Eqs.  (13)  to  (15)  to 
calculate  A[Hb02],  A[Hb],  and  A[Hbtotal],  A[Hblotal]  re- 
mained  constant  during  the  oxygenation  and  deoxygenation 
cycles  as  expected.  This  demonstrates  that  the  values  of 
Pmoi=  1*103  and  jSHb=  0.9035  are  correct  and  necessary  to 
compensate  the  differences  in  DPFs  caused  by  the  two  differ¬ 
ent  wavelengths. 

3,4  Tissue  Phantom  Solution  Model 

In  order  to  study  the  relationship  between  p02  and  A[Hb02] 
in  regular  tissues,  we  conducted  a  tissue-simulating  phantom 
study  by  using  the  liquid  solution  similar  to  that  mentioned 
above.  In  normal  tissues,  there  are  several  steps  of  oxygen 
transport  from  the  blood  to  tissue  cells.23  In  the  tissue- 
simulating  phantom,  blowing  oxygen  gas  represents  oxygen¬ 
ation  process  of  blood  in  the  lungs,  and  blowing  nitrogen  gas 
simulates  deoxygenation  process  of  blood  in  the  tissues.  The 


%  2  Simultaneous  dynamic  changes  of  A[Hb02],  A[Hb],  and 
A[Hb]to{3|  in  the  phantom  solution  measured  using  NJRS.  The  gray 
solid  curve  is  for  A[Hb]tolaI  without  using  /?Hbo2  and  finb  values.  Oxy¬ 
gen  consumption  by  yeast  produced  deoxygenated  blood  and  blow¬ 
ing  oxygen  restored  oxygenation.  During  the  oxy-  and  deoxygenation 
process,  A[Hb],otaf  is  supposed  to  be  a  constant.  However,  as  we  can 
see  here,  &[Hb]lot3j  calculated  without  /?Hb03  arjd  0m  values  shows 
the  fluctuation  during  the  oxy-  and  deoxygenation  while  A[Hb]tota| 
calculated  using  $Hb02  ar|d  $Hb  values  shows  the  veracity  of  these 
modified  algorithms. 


differences  between  the  tissue-simulating  phantom  and  real 
tissues  are  that  there  is  no  capillary  membrane  in  the  phan¬ 
tom,  and  that  the  phantom  is  more  homogeneous  than  real 
tissues.  Capillary  membranes  have  high  permeability  of  oxy¬ 
gen,  so  oxygen  transport  from  blood  to  tissues  crossing  the 
capillary  membranes  occurs  straightforwardly.  Furthermore, 
normal  tissues  are  well  vascularized,  and  the  MR  techniques 
are  more  sensitive  toward  measuring  small  vessels  and  vascu¬ 
lar  bed  of  the  tissue.24  Therefore,  vasculature  of  normal  tis¬ 
sues  has  been  simulated  by  a  turbid  solution  mixed  with  blood 
as  a  simplified  laboratory  model  in  NIRS  measurements  for 
oxygen  transport  from  blood  to  normal  tissues.10,18,22 

The  experimental  setup  shown  in  Figure  3  was  made  to 
simulate  tumor  oxygenation/deoxygenation.  Oxygen  needle 
electrodes,  a  pH  electrode,  and  a  thermocouple  probe  (model 
2001,  Sentron,  Inc.,  Gig  Harbor,  WA)  were  placed  in  the  so¬ 
lution,  and  the  gas  tube  for  delivery  of  N2  or  air  was  placed 
opposite  the  NIRS  probes  to  minimize  any  liquid  movement 
effects.  Source  and  detector  probes  for  the  NIRS  were  placed 
in  reflection  geometry  with  a  direct  separation  of  3  cm.  The 
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Fig.  3  Experimental  setup  for  phantom  study  using  1%  Intralipid  in 
saline  buffer.  NIRS  probes  were  placed  in  reflectance  mode,  while  the 
gas  bubbler  was  placed  opposite  to  minimize  liquid  movement  ef¬ 
fects.  After  adding  2  ml  of  rabbit  blood  to  a  200  ml  solution,  nitrogen 
gas  and  air  were  introduced  to  deoxygenate  and  oxygenate  the  solu¬ 
tion,  respectively. 


journal  of  Biomedical  Optics  *  January  2003  *  Vol.  8  No.  1  57 


solution  was  stirred  constantly  to  maintain  homogeneity  by  a 
magnetic  stirrer  at  ~37  °C  Fresh  whole  rabbit  blood  (2  mL) 
was  added  to  the  200  mL  solution  before  baseline  measure¬ 
ment.  Nitrogen  gas  and  air  were  used  to  deoxygenate  and 
oxygenate  the  solution,  respectively. 

3.5  MRI  Instrumentation  and  Procedure 

To  support  the  findings  obtained  from  the  p02  electrode  mea¬ 
surements  and  MRS,  we  conducted  MRI  experiments  using 
an  Omega  CSI  4.7  T  40  cm  system  with  actively  shielded 
gradients.  A  homebuilt  tunable  1H/19F  single  turn  solenoid 
coil  was  placed  around  the  tumor.  45  fib  hexafluorobenzene 
(HFB;  Lancaster,  Gainesville,  FL)  was  administered  directly 
into  the  tumor  using  a  Hamilton  syringe  (Reno,  NV)  with  a 
custom-made  fine  sharp  (32  gauge)  needle  and  HFB  was  de¬ 
liberated  dispersal  along  several  tracks  to  interrogate  both 
central  and  peripheral  tumor  regions,  as  described  in  detail 
previously.5  HFB  is  ideal  for  imaging  p02  because  it  has  a 
single  resonance  and  its  relaxation  rate  varies  linearly  with 
oxygen  concentration.  ]H  images  were  acquired  for  anatomi¬ 
cal  reference  using  a  traditional  3-D  spin-echo  pulse  se¬ 
quence.  Conventional  19F  MR  images  were  taken  to  show  the 
3-D  distribution  of  the  HFB  in  the  tumor.  19F  MR  images 
were  directly  overlaid  over  !H  images  to  show  the  position  of 
the  HFB  in  that  slice. 

Tumor  oxygenation  was  assessed  using  fluorocarbon  relax- 
ometry  using  echo  planar  imaging  for  dynamic  oxygen  map¬ 
ping  (FREDOM)  based  on  l9F  pulse  burst  saturation  recovery 
(PBSR)  echo  planar  imaging  (EPI)  of  HFB.25  The  PBSR 
preparation  pulse  sequence  consists  of  a  series  of  20  nonspa- 
tially  selective  saturating  90  deg  pulses  with  20  ms  spacing  to 
saturate  the  19F  nuclei.  Following  a  variable  delay  time  %  a 
single  spin-echo  EPI  sequence  with  blipped  phase  encoding 
was  applied.26  Fourteen  32X32  PBSR-EPI  images,  with  r 
ranging  from  200  ms  to  90  sec  and  a  field  of  view  (FOV)  of 
40X40  mm,  were  acquired  in  8  min  using  the  alternated  re¬ 
laxation  delays  with  variable  acquisitions  to  reduce  clearance 
effects  (ARDVARC)  acquisition  protocol.25  An  Rl(=  1/71 ) 
map  was  obtained  by  fitting  the  signal  intensity  of  each  voxel 
of  the  14  images  to  a  three-parameter  relaxation  model  by  the 
Levenberg-Marquardt  least-squares  algorithm: 

y*(U)=A(iJH  1  -(1  +  W)zxp(  —  Rl(iJ)-Tn)} 

(24) 


(»  =  1,2,..  • 

••■,14) 

•••,32), 

where  yn(ij )  is  the  measured  signal  intensity  corresponding 
to  delay  time  rn  (the  w’th  images)  for  voxel  (i,j\  A  (ij)  is  the 
fully  relaxed  signal  intensity  amplitude  of  voxel  (i,  J)9  W  is  a 
dimensionless  scaling  factor  allowing  for  imperfect  signal 
conversion,  Rl(iJ)  is  the  relaxation  rate  of  voxel  (i,  j)  in 
units  of  sec-1,  and  A,  W}  and  Rl  are  the  three  fit  parameters 
for  each  of  the  32X32  voxels.  Finally,  the  p02  maps  were 
generated  by  applying  the  calibration  curve,  p02  (mm  Hg) 
=  [i?  1  (s“ !)— 0.0835]/0.001876  at  37°C,  to  the  Rl  maps.25 
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Fig*  4  Simultaneous  dynamic  changes  of  A[Hb02]  and  p02  in 
R3327-H!  rat  prostate  tumors  using  N1RS  and  p02  needle  electrode, 
(a)  A  small  tumor  {1.5  cm3)  showed  a  rapid  p02  response  {case  1)# 
whereas  (b)  a  bigger  tumor  (3.1  cm3)  showed  a  slower  p02  response 
(case  2),  (c)  In  a  third  tumor  (1.6  cm3)  where  regional  baseline  p02 
was  <5  mm  Hg,  there  was  no  p02  response  (case  3).  The  unit  of 
A[Hb02]  is  mM/DPF,  where  DPF  is  equal  to  the  optical  path  length 
divided  by  the  source-detector  separation.  Dotted  vertical  line  marks 
the  time  when  the  gas  was  changed. 

4  Results 

4.1  Tumor  Study  Results 

We  have  measured  relative  changes  of  [Hb02],  [Hb], 
[Hbtotat],  and  tumor  tissue  p02  (electrode)  from  eight  Dun¬ 
ning  prostate  R3327-HI  tumors,  and  Figure  4  shows  three 
representative  data  sets.  Figure  4(a)  shows  the  temporal  pro¬ 
files  of  A[Hb02]  and  p02  in  a  small  Dunning  prostate 
R3327-HI  tumor  (1.5  cm3)  measured  simultaneously  with 
NIRS  and  the  p02  needle  electrode  during  respiratory  chal¬ 
lenge.  After  a  switch  from  air  to  carbogen,  A[Hb02]  in¬ 
creased  rapidly,  along  with  tumor  tissue  p02 .  Figure  4(b)  was 
obtained  from  a  large  tumor  (3.1  cm3):  the  electrode  readings 
showed  a  slower  p02  response,  whereas  the  NIRS  response 
was  biphasic,  which  has  been  a  commonly  observed  dynamic 
feature.11  In  a  third  tumor  (1.6  cm3),  NIRS  behaved  as  before, 
but  p02  did  not  change  [Figure  4(c)]. 


Interplay  of  Tumor  Vascular  Oxygenation  . . . 


(a)  Air  Carbogen 


Fig.  5  Dynamic  changes  of  A[Hb02]  and  p02  in  two  R3327-AT1  rat 
prostate  tumors  measured  sequentially  using  NtRS  and  19F  MR  p02 
mapping.  The  solid  curves  represent  A[Hb02],  and  the  solid  lines 
with  solid  circles  represent  mean  p02±$E  (standard  error)  of  21  (a) 
and  45  (b)  voxels  of  the  respective  tumor.  Dashed  lines  with  open 
symbols  are  4  representative  voxels  for  each  case.  After  a  gas 
switched  from  air  to  carbogen,  the  mean  p02  values  of  both  tumors 
increased,  but  individual  voxels  showed  quite  different  responses,  in¬ 
dicating  high  heterogeneity  in  the  tumors.  The  tumor  sizes  were  3.2 
cm3  and  2.7  cm3  for  (a)  and  (b),  respectively. 


In  four  tumors  from  a  separate  subline  (Dunning  prostate 
R3327-AT1),  NIRS  and  l9F  MRI  were  taken  sequentially  with 
carbogen  challenge,  and  two  representative  data  sets  are 
shown  in  Figure  5.  NIRS  response  showed  vascular  oxygen¬ 
ation  changes  as  before,  and  FREDOM  revealed  the  distinct 
heterogeneity  of  the  tumor  tissue  response.  Initial  pC>2  was  in 
the  range  of  1  to  75  mm  Hg,  and  carbogen  challenge  produced 
p02  values  in  the  range  of  6  to  350  mm  Hg.  Representative 
voxels  are  shown  in  each  figure  by  dashed  lines  with  open 
symbols.  In  addition,  mean  p02  values  were  calculated  by 
averaging  all  available  p02  readings  over  21  and  45  voxels 
for  the  two  respective  tumors.  We  usually  obtain  p02  tempo¬ 
ral  profiles  from  individual  voxels  among  200  to  400  voxels 
in  a  tumor  during  the  entire  intervention  period.  The  p02 
readings  presented  here  were  picked  to  show  heterogeneity  of 
the  tumor.  In  Figure  5(a),  the  closest  distance  between  the  two 
voxels  is  L25  mm  (between  0  and  □),  and  the  furthest  dis¬ 
tance  is  7.6  mm  (between  X  and  A).  In  Figure  5(b),  the 
closest  distance  is  3.6  mm  (between  X  and  A)  and  the  fur¬ 
thest  distance  is  16  mm  (between  X  and  □).  These  indeed 
showed  that  tumor  p02  responses  to  carbogen  intervention 
could  be  quite  different  at  different  locations.  Notice  that  Fig¬ 
ure  5(a)  showed  spikes  of  A[Hb02]  during  the  measurement. 
We  expect  this  to  be  caused  by  sudden  changes  in  rat  respi¬ 
ratory  circulation  or  motion,  rather  than  resulting  from  simple 
instrumental  noise.  It  is  also  seen  that  mean  p02  values  have 
displayed  a  consistent  increase  when  A[Hb02]  showed 
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Fig.  6  Simultaneous  dynamic  changes  of  A[Hb02],  A[Hb]fola|,  and 
p02  in  the  phantom  solution  measured  using  NIRS  and  p02  needle 
electrode.  The  dark  solid  curve  is  for  A[Hb02],  the  lighter  solid  fine  is 
for  A[Hb]tola(#  and  the  solid  circles  show  p02  values  in  the  phantom 
solution.  After  ~3  mm  baseline,  the  bubbling  gas  was  changed  from 
air  to  nitrogen  to  deoxygenate  the  solution  and  then  switched  back  to 
air  to  reoxygenate  the  solution.  The  unit  of  A[Hb02]  is  mM/DPF. 


spikes,  suggesting  that  such  spikes  may  result  from  changes  in 
rat  physiological  conditions. 

4.2  Tissue  Phantom  Study  Results 

Figure  6  shows  a  temporal  profile  for  A[Hb02]  and  p02  mea¬ 
sured  from  the  tissue  phantom  during  a  cycle  of  gas  change 
from  air  to  nitrogen  and  back.  The  first  three  minutes  were 
measured  as  a  baseline  after  adding  2  ml  blood.  Bubbling 
nitrogen  deoxygenated  the  solution  and  caused  the  p02  values 
to  fall;  A[Hb02]  declined  accordingly  with  a  small  time  lag. 
After  the  bubbling  gas  was  switched  from  nitrogen  to  air,  both 
A[Hb02]  and  p02  started  to  increase  simultaneously,  but  the 
recovery  time  of  A[Hb02]  to  the  baseline  was  faster  than  that 
of  p02 ,  The  small  time  lag  between  the  changes  of  A[Hb02] 
and  p02  is  probably  due  to  the  allosteric  interactions  between 
hemoglobin  and  oxygen  molecules.  According  to  the  hemo¬ 
globin  oxygen-dissociation  curve,19127  oxyhemoglobin  starts 
to  lose  oxygen  significantly  when  p02  falls  below  70  mm  Hg 
at  standard  conditions  (pH  =7.4,  pCO2=40mmHg,  and 
temperature=37°C).  The  same  principle  can  explain  why 
A[Hb02]  has  a  faster  recovery  than  that  of  p02.  Figure  6 
shows  that  A[Hb02]  is  already  saturated  when  p02  is  at  50 
mm  Hg,  while  the  solution  was  still  being  oxygenated.  This 
may  be  due  to  low  pC02  in  the  solution  where  this  can  shift 
the  oxyhemoglobin  dissociation  curve  to  the  left,  causing  oxy¬ 
hemoglobin  to  be  saturated  at  lower  p02.  Importantly, 
A  [Hb]  totai  remained  unchanged,  as  expected,  during  a  cycle 
of  deoxygenation  and  oxygenation. 

43  Correlation  between  p02  and  Normalized 
A[Hb02] 

For  Tissue  Phantoms .  Figure  7(a)  replots  the  data  given  in 
Figure  6,  showing  the  relationship  between  normalized 
A[Hb02]  and  p02  measured  from  the  tissue  phantom  during 
the  oxygenation  (air  blowing)  period  after  the  nitrogen  blow¬ 
ing.  Open  circles  are  the  measured  data,  .and  the  solid  line  is 
the  fitted  curve  using  Eq.  (23).  The  error  bars  for  the  data 
were  not  shown  here  since  they  are  smaller  than  the  symbols 
of  the  data  points.  For  the  curve  fitting  procedure,  we  used  a 
nonlinear  curve-fitting  routine  provided  through  Kaleida- 
Graph  (Synergy  software,  Reading.  PA),  The  fitted  parameters 
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Tumor  tissue  p02  (mmHg) 

Fig.  7  Changes  of  tissue  p02  with  normalized  changes  of  oxygenated 
hemoglobin  (a)  in  the  phantom  solution  using  the  N1RS  and  p02 
needle  electrode  and  (b)  in  tumors  measured  with  NfRS,  p02  needle 
electrode,  and  19F  MR  p02  mapping.  In  (a),  the  open  circles  are  mea¬ 
sured  data  and  the  solid  line  is  the  fitted  curve  using  Eq.  (21).  This 
shows  that  Eq.  (21)  works  well  in  a  homogeneous  system.  In  (b),  all 
the  tumor  data  are  shown  indicating  that  tumors  are  highly  heteroge¬ 
neous  for  p02  response  to  carbogen  inhalation.  Open  symbols  show 
local  p02  changes  (from  Figure  4)  and  solid  symbols  show  the  mean 
p02  changes  (from  Figure  5)  during  gas  intervention.  To  estimate  glo¬ 
bal  s02  in  tumors  during  respiratory  challenges,  we  applied  Eq.  (23) 
to  Figure  5(a),  indicating  s02  changes  during  carbogen  inhalation 
when  compared  via  tumor  pD2 . 


are  n=  1.9,  P50=  15.2  mm  Hg,  lsO2’]\msc-0%t  and 

[s02]max”99%  with  R  —  0,997  and  minimized  chi-square. 
The  fitted  values  of  [s02]base  and  [s02]max  are  in  good  agree¬ 
ment  with  the  expected  values,  since  the  corresponding  p02 
values  are  0  and  160  mm  Hg,  respectively.  This  agreement 
validates  Eq.  (23)  and  further  indicates  that  we  can  measure 
approximate  s02  values  during  the  gas  interventions  in  a  ho¬ 
mogeneous  system  by  fitting  the  experimental  data  using  Eq. 
(23)  even  though  we  do  not  measure  absolute  [Hb02].  The 
Hill  coefficient  (n)  and  p02  value  at  50%  of  s02  (P50)  are 
smaller  than  the  values  from  a  standard  oxyhemoglobin  satu¬ 
ration  curve,  probably  due  to  the  shift  of  the  oxyhemoglobin 
dissociation  curve. 

For  Tumor  Study.  Figure  7(b)  replots  the  data  given  in  Figures 
4  and  5,  showing  a  direct  relationship  between  the  normalized 
A[Hb02]  and  tissue  p02  in  the  tumors.  NIR3  results  tended 
to  be  similar  for  several  tumors,  and  p02  electrode  measure¬ 
ments  showed  considerable  variation  even  in  the  same  tumor 
type,  suggesting  distinct  tumor  heterogeneity.  This  was  sub¬ 
stantiated  by  the  I9F  MR  p02  mappings  (Figure  5);  indeed,  in 
some  cases,  p02  values  did  not  change  with  respiratory  chal¬ 
lenge.  especially  when  baseline  p02  values  were  lower  than 
10  mm  Hg. 

Equation  (23)  can  be  used  to  estimate  values  of  [s02]base 
and  [s02]max  for  the  tissue-simulating  phantom  (a  homoge¬ 
neous  system).  However,  the  relationship  fails  for  heteroge- 
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neous  systems  such  as  tumors.  The  NIRS  measurements  in¬ 
terrogate  a  large  volume  of  tumor  tissue,  giving  a  global  value 
of  normalized  A[Hb02],  whereas  the  p02  readings  are  local 
near  the  tip  of  the  needle  electrode.  However,  to  estimate 
mean  values  of  [s02]basc  and  [s02]max>  R  is  reasonable  to 
compare  the  global  normalized  A[Hb02]  with  global  tissue 
p02 ,  which  can  be  obtained  by  summing  up  all  local  p02 
readings  at  different  pixels  measured  from  the  19F  MRI  map¬ 
ping,  as  done  in  Sec,  4,1  and  shown  by  solid  lines  in  Figure  5. 
The  data  shown  in  Figure  7(b)  with  solid  symbols  are  the 
global  mean  p02  values  calculated  from  the  corresponding 
MRI  data.  The  solid  fitting  curve  shown  in  Figure  7(a)  is 
obtained  from  the  mean  p02  data  given  in  Figure  5(a),  In  this 
case,  the  fitting  parameters  are  P50,  [s02]base,  and  [s02]max 
with  a  fixed  Hill  coefficient  n  to  be  the  same  as  that  under 
standard  conditions.  The  best  fitting  curve  of  Eq,  (23)  is 
shown  in  Figure  7(b),  having  P5O=20,6±4,l  mm  Hg, 
[s02]base=37±  13%,  and  [s02]max— 100%  with  i?  =  0.985 
and  goodness-of-fit  %2— 0.031.  Estimated  errors  for  P50  and 
[s02]base  arc  not  insignificant  and  a  better  fit  could  be  found 
by  measuring  p02  with  better  temporal  resolution. 

5  Discussion  and  Conclusion 

Tumor  oxygenation  involves  a  complex  interplay  of  multiple 
compartments  and  parameters:  blood  flow,  blood  volume, 
blood  vessel  structure,  and  oxygen  consumption.  NIRS  pro¬ 
vides  a  global  noninvasive  estimate  of  average  vascular  oxy¬ 
genation  encompassing  arterial,  venous,  and  capillary  com¬ 
partments,  In  agreement  with  our  previous  observations,11  the 
A[Hb02]  response  is  often  biphasie,  which  we  believe  repre¬ 
sents  rapid  elevation  of  arterial  oxygenation,  followed  by 
more  sluggish  capillary  components. 

Comparison  with  simultaneous  electrode  measurements  in¬ 
deed  revealed  that  tumors  are  heterogeneous.  Like  NIRS  mea¬ 
surements,  p02  electrodes  provide  rapid  assessment  of  p02 
facilitating  real-time  observation  of  dynamic  changes.  In  Fig¬ 
ure  4(a),  p02  starts  at  a  baseline  value  —15  mm  Hg  and 
increases  rapidly  in  response  to  respiratory  challenge  with 
carbogen.  Indeed,  the  rate  approaches  that  of  the  vascular 
compartments.  In  a  second  tumor  [Figure  4(b)],  where  the 
interrogated  location  showed  a  slightly  lower  p02 ,  the  tissue 
response  was  more  sluggish.  For  a  third  HI  tumor,  local  base¬ 
line  p02  was  found  to  be  <5  mm  Hg,  and  this  did  not  change 
with  carbogen  inhalation  despite  the  response  observed  by 
NIRS.  This  suggests  a  danger  of  comparing  a  global  vascular 
measurement  with  regional  tumor  p02,  since  tumors  are 
known  to  be  highly  heterogeneous.  This  also  demonstrates  an 
essential  need  for  NIR  imaging  of  tumors  to  provide  regional 
tumor  vascular  oxygenation  details, 

FREDOM  measurements  in  Figure  5  revealed  the  hetero¬ 
geneity  in  baseline  oxygenation  within  individual  tumors  of 
this  second  tumor  subline  as  also  reported  previously,25  Base¬ 
line  p02  ranged  from  1  to  75  mm  Hg,  and  response  to  carbo¬ 
gen  was  variable  in  terms  of  rate  and  extent,  as  also  seen  for 
the  HI  subline  using  electrodes  (Figure  4).  As  with  the  elec¬ 
trodes,  the  better  oxygenated  tumor  regions  showed  a  faster 
and  greater  response  to  carbogen  inhalation.  The  oxygen  elec¬ 
trode  measurements  in  Figure  4  showed  a  maximum  p02  of 
around  45  mm  Hg,  though  we  have  observed  values  as  high  as 
95  mm  Hg  using  oxygen  needle  electrode.  Observations  using 
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the  fluorescence-based  OxyLite™  fiber-optic  devices  for  mea¬ 
suring  HI  tumor  reached  the  maximum  detectable  p02  of  100 
mm  Hg  during  carbogen  inhalation.5  FREDOM  has  shown 
values  of  less  than  5  mm  Hg  and  greater  than  160  mm  Hg 
under  air  breathing  conditions,  and  reaching  350  mm  Hg  in 
HI  tumors  while  breathing  carbogen,5  Each  method  indicates 
that  tumors  are  highly  heterogeneous,  but  it  has  been  shown 
that  there  can  be  a  positive  linear  relationship  between  base¬ 
line  p02  and  maximum  p02  during  carbogen  inhalation  in  the 
Dunning  prostate  ATI  tumor  line,8 

The  phantom  measurements  indicate  and  validate  the  reli¬ 
ability  of  the  NIRS  technique  and  also  prove  that  normalized 
A[Hb02]  is  closely  related  to  the  normalized  hemoglobin- 
oxygen  dissociation  curve.  The  phantom  data  confirmed  that 
we  can  obtain  absolute  s02  values  in  a  homogeneous  system 
by  measuring  both  A[Hb02]  and  p02.  We  could  estimate 
mean  s02  values  of  the  tumor  under  intervention  using  global 
A[Hb02]  and  averaged  p02  readings,  and  the  fitting  errors 
are  expected  to  be  improved  by  having  more  data  points. 
Measuring  regional  tumor  vascular  oxygenation  by  NIR  im¬ 
aging  of  tumors  should  allow  us  to  correlate  local  A[Hb02] 
and  p02  and  to  understand  the  oxygen  transport  process  from 
tumor  vasculature  to  tumor  tissue,  and  this  is  the  direction  of 
our  future  work. 

Both  NIRS  and  electrodes  offer  essentially  real-time  mea¬ 
surement  of  changes  in  oxygenation,  which  can  be  rapid  (Fig¬ 
ure  4),  Indeed,  the  inflow  kinetics  of  vascular  02  detected  by 
NIRS  are  similar  to  those  previously  reported  in  the  HI  tumor 
line  following  a  bolus  of  the  paramagnetic  contrast  agents 
Gd-DTPA.28  FREDOM  has  lower  temporal  resolution,  but  re¬ 
veals  the  tumor  heterogeneity  and  differential  response  of  re¬ 
gions  exhibiting  diverse  baseline  p02 .  The  results  here  cor¬ 
respond  closely  with  more  extensive  observation.5,8*25  While 
FREDOM  currently  requires  6.5  min  per  p02  map,  we  have 
previously  demonstrated  an  alternative  data  acquisition  proto¬ 
col  achieving  1  s  time  resolution  in  a  perfused  heart,  albeit 
providing  less  precision  in  measurements  and  only  a  global 
determination  29  Such  an  approach  could  allow  us  to  measure 
global  A[Hb02]  and  global  p02  simultaneously  with  a  high 
temporal  resolution,  understand  the  relationship  between  glo¬ 
bal  A[Hb02]  and  global  p02,  and  obtain  absolute  values  of 
s02  of  the  tumors  as  tumors  grow. 

In  conclusion,  we  have  refined  the  algorithms  for  calculat¬ 
ing  [Hb],  [Hb02],  and  [Hbtotai]  and  measured  relative 
[Hb02]  changes  in  tumor  vasculature  and  tumor  tissue  p02 
under  carbogen  intervention  using  NIRS  and  a  needle  type 
p02  electrode.  The  p02  data  were  also  supported  by  the  t9F 
MR  p02  mapping.  We  have  also  developed  an  algorithm  to 
estimate  s02  values  in  the  tumor  during  respiratory  interven¬ 
tions.  The  NIRS  data  showed  significant  changes  in  vascular 
oxygenation  accompanying  respiratory  interventions,  and 
changes  in  tumor  vascular  oxygenation  preceded  tumor  tissue 
p02 .  Oxygen  electrode  measurements  and  19F  MR  p02  map¬ 
ping  results  proved  that  tumors  are  highly  heterogeneous.  The 
phantom  data  confirmed  that  normalized  [Hb02]  data  to¬ 
gether  with  p02  measurements  can  be  used  to  estimate  abso¬ 
lute  s02  values  in  a  homogeneous  system.  For  a  highly  het¬ 
erogeneous  medium,  such  as  tumors,  local  comparison 
between  the  [Hb02]  and  p02  value  is  desired  and  required  in 
order  to  reveal  the  process  of  oxygen  delivery  from  the  tumor 


vascular  bed  to  the  tumor  tissues.  Therefore,  this  study  not 
only  demonstrates  that  the  NIRS  technology  can  provide  an 
efficient,  real-time,  noninvasive  approach  to  monitoring  tumor 
physiology  and  is  complementary  to  other  techniques,  but 
also  emphasizes  the  need  to  develop  an  imaging  technique  to 
study  spatial  heterogeneity  of  tumor  vasculature  under  oxygen 
or  other  therapeutic  interventions. 
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